NASA Contractor Report 191087 


/A 

/ 5 <& 3/7 

jut? 


Unsteady Blade Pressures on a Propfan: Predicted 
and Measured Compressibility Effects 


M. Nallasamy 
Sverdrup Technology , Inc. 
Lewis Research Center Group 
Brook Park, Ohio 


July 1992 


( NASA-C R- 19108 7 ) UNSTEADY BLADE N93-23431 

PRESSURES ON A PROPFAN: PREDICTED 
AND MEASURED COMPRESSIBILITY 

EFFECTS Final Report (Sverdrup Unclas 

Technology) 16 p 


G3/02 0156317 


Prepared for 

Lewis Research Center 

Under Contract NAS3-25266 


NASA 

National Aeronautics and 
Space Administration 


UNSTEADY BLADE PRESSURES ON A PROPFAN: 
PREDICTED AND MEASURED COMPRESSIBILITY EFFECTS 


M. Nallasamy 
Sverdrup Technology, Inc. 
NASA Lewis Research Center Group 
Cleveland, Ohio 


Abstract 

The effect of compressibility on unsteady blade 
pressures is studied by solving the three-dimensional 
Euler equations. The operation of the eight-bladed 
SR7L propfan at 4.75° angle of attack was consid- 
ered. Euler solutions were obtained for three Mach 
numbers, 0.6, 0.7 and 0.8 and the predicted blade 
pressure waveforms were compared with flight data. 
The comparisons show that in general, the effect of 
Mach number on pressure waveforms are correctly 
predicted. The change in pressure waveforms are min- 
imal when the Mach number is increased from 0.6 to 
0.7. Increasing the Mach number from 0.7 to 0.8 pro- 
duces significant changes in predicted pressure lev- 
els. The predicted amplitudes, however, differ from 
measurements at some transducer locations. Also the 
predicted appearance of a shock in the highly loaded 
portion of the blade revolution is not indicated by 
the measurements. At all the three Mach numbers, 
the measured (installed propfan) pressure waveforms 
show a relative phase lag compared to the computed 
(propfan alone) waveforms due to installation effects. 
Measured waveforms in the blade tip region show 
nonlinear variations which are not captured by the 
present numerical procedure. 

Introduction 

The flow unsteadiness affects critically the prop- 
fan performance and near field noise levels. The un- 
steadiness may be the r^ult of operation of the prop- 
fan ^a^ an angle relative to The mean flow direction 
or of distorted inflow caused by installation effects. 
Wind tunnel and flight tests were conducted on a 2.74 
m (9 ft) diameter large-scale propfan to further un- 
derstand the effects of flow unsteadiness 1-3 . 

Unsteady blade pressure measurements of the 
large-scale propfan were first made in a transonic 
wind tunnel with angular inflow and (cylinder) wake 
inflow 1 . It was found that at takeoff conditions, high 
power cases resulted in the formation of a leading edge 
vortex on the blade. Then in the Propfan Test Assess- 
ment (PTA) program, flight tests were conducted to 
investigate the effect of inflow angle on the near-field 
noise level 2 . A sensitivity of about ldB per degree of 


inflow angle change was found in flight measurements. 
Finally in the PTA follow-on flight tests, detailed un- 
steady blade pressure measurements 3 were made on 
a specially designed instrumented blade. In this test, 
the blade suction surface had 20 pressure transducers 
distributed over three radial stations ( r/R = 0.68, 
0.86 and 0.95, where r is the radial distance and R is 
the blade tip radius) while the pressure surface had 
10 pressure transducers distributed over two radial 
stations ( r/R = 0.68 and 0.95). A nacelle tilt ar- 
rangement was employed to vary the inflow angle to 
the propfan. 

Three sets of data were taken providing an 
unique database of detailed unsteady blade pressures 
in flight: 1. Low altitude (580 m), low speed (Mach 
number ~ 0.3) 2. High altitude (10,500 m), high 
speed (Mach number ~ 0.8) 3. A compressibility se- 
ries in which the flight Mach number ( M ) was varied 
(Mach number ~ 0.6, 0.7 and 0.8), keeping the ad- 
vance ratio and power coefficient (blade angle) con- 
stant. In each set, three nacelle tilt angles, -3, -1 (tilt 
down) and 2° (tilt up), were considered giving an ef- 
fective inflow angle variation of 5°. Efforts have been 
made to understand the unsteady flow features exhib- 
ited by these data and also to validate the prediction 
of three dimensional Euler analysis methods against 
the data. 

The unsteady flow features of a propfan at 
takeoff (low speed, low altitude) were examined by 
Nallasamy 4 by solving the three-dimensional Euler 
equations. It was found that the measured (installed 
propfan) blade pressure waveforms showed a relative 
phase lag compared to the computed (propfan alone) 
waveforms at the outboard ( r/R = 0.95) station on 
both pressure and suction surfaces. However, on the 
pressure surface the magnitudes were found to be in 
good agreement with flight data at all inflow angles 
studied. On the suction surface, in addition to the rel- 
ative phase lag. the measurements showed distortion 
of the waveforms. The extent of distortion increased 
with inflow angle. 

Nallasamy and Groeneweg 5 studied the unsteady 
flow features at cruise operating conditions. They 
showed that at inflow angles of 1.6 and 4.6 degrees 
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passage shocks extending from suction to pressure 
surface formed and dissolved during each revolution 
of the blade. The computed unsteady blade pressures 
for 1.6° were compared with wind tunnel data 6 . The 
comparisons showed good agreement of the predicted 
blade pressure waveforms with data at most of the 
transducer locations considered. 

The effect of compressibility on steady blade 
pressures was examined 7 . It was found that with in- 
creasing Mach number, more significant changes oc- 
curred on the suction side than on the pressure side. 
The evidence of a compression wave started to de- 
velop when the Mach number reached 0.7. At a free 
stream Mach number of 0.78, the compression wave 
fully developed into a trailing edge shock. 

In the present paper, the effects of compressibil- 
ity on unsteady blade pressures of a propfan are stud- 
ied by solving the three dimensional Euler equations 
for angular inflow. The solutions are obtained for 
three Mach numbers, 0.6, 0.7 and 0.8 of the PTA 
follow-on compressibility test series, keeping advance 
ratio and blade angle constant. First, the predicted 
blade pressure waveforms for M = 0.8 are compared 
directly with PTA follow-on flight data to assess the 
predictions at the high inflow angle, high speed case. 
Then the waveforms for three Mach numbers, 0.6, 0.7 
and 0.8, are compared with measured waveforms to 
evaluate the ability of the solution procedure in pre- 
dicting the Mach number effects observed in flight 
tests. No effort is made here to quantify the effect of 
compressibility on propfan performance. 

Numerical Solution of Three-Dimensional 
Euler Equations 

The unsteady three-dimensional Euler equations 
governing the inviscid flow through a propfan are 
solved employing a solution procedure developed by 
Whitfield et al. 8,9 . In this procedure the Euler equa- 
tions in conservative differential form are transformed 
from a Cartesian reference frame to a body fitted 
curvilinear reference frame. The transformed equa- 
tions are then discretized employing a finite volume 
technique. A Lower-Upper (LU) implicit numerical 
scheme is used to solve the discretized equations. An 
approximate Riemann solver is used for block inter- 
face definitions. The flowfield is represented by a 
multi-block composite grid to limit the core memory 
requirements. 

Flow Configuration and Computational Grid 

The configuration considered is that of the eight 
bladed SR7L propfan of the flight test 3 . The direction 
of rotation of the propfan is clockwise, looking down- 
stream. The azimuth angle, $, is measured from ver- 


tical (top-dead-center) for aircraft installation as in 
the presentation of flight data in Ref.3 and increases 
in the direction of rotation (Fig. 0). The grid em- 
ployed is an H-grid with 107 x 41 x 25 (axial by radial 
by circumferential) nodal points in each passage. The 
blade design hot coordinates are used to generate the 
grid. Each passage is divided into three blocks with 
107 x 41 x 9 grid points in each block. Thus 24 blocks 
of grid were used to describe the entire flowfield. Each 
blade surface is represented by 49 x 27 (axial by ra- 
dial) grid points with higher resolution near the lead- 
ing and trailing edges, the hub and the tip. 

Results and Discussion 

The unsteady three-dimensional Euler solutions 
were obtained for the following three flight test cases 
of Ref. 3 : (i) run = 144, id = lc and Mach number, 
M — 0.603, (ii) run = 142, id = 2c, M — 0.701 and 
(iii) run = 140, id = 3c and M — 0.803. These three 
runs had a nacelle tilt of +2°. The effective inflow 
angle is dependent not only on the nacelle tilt angle 
but also on the airplane angle of attack and the up- 
wash angle at the propfan (Fig. 0). For the above test 
runs the average value of the airplane angle of attack 
was 1.75°. The estimated upwash 5 at the propfan was 
1.0°. Thus the effective inflow angle to the propfan 
was 4.75°. The Euler solutions were obtained with an 
advance ratio of 3.17, blade setting angle of 57° and 
inflow angle of 4.75° for all three Mach numbers. 

The Euler solutions are obtained from an impulse 
start for three complete revolutions of the blade, to 
obtain a reasonably accurate solution. The results 
of the third revolution are analyzed and presented 
here. The predicted total power coefficient (for eight 
blades) are shown in Table I along with the measured 
ones. In the experiments, the blade angle for M — 0.6 
is different from the other two. In the predictions, 
though the blade angle (57°) is maintained the same 
for all Mach numbers, the total power coefficient is 
under predicted by about 6% for M = 0.6. The total 
power coefficients for M = 0.7 and 0.8 are in rea- 
sonable agreement with data. The power per blade 
variation was found to show the expected sinusoidal 
variation due to inflow angle, in all three cases. The 
amplitude of the stabilized power coefficient during 
the third revolution varies -f 109% and -98% about the 
mean for M - 0.6, +120% and -102% for M = 0.7 
and +105% and -104% for M — 0.8. 

First, the predicted and measured blade pressure 
waveforms are compared for the high Mach number 
case, M = 0.8. Figure 1 shows the blade pressure 
waveforms at the transducer locations ( xfc , where x is 
the axial distance and c is the blade chord) on the suc- 
tion surface at the inboard radial station r/R = 0.68. 
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The predicted amplitudes agree fairly well with flight 
data. The measured waveforms show a relative phase 
lag compared to the predicted ones. This is in con- 
trast to the takeoff case studied 4 where such a phase 
lag was observed only at the outboard radial station. 
The observed phase lag seems to be due to installation 
effect, that is, the presence of the wing in flight tests 
as compared to the propfan alone configuration of 
the computation. Heidelberg and Woodward 10 noted 
similar phase variations in their model tests in the 
wind tunnel with and without wing installation. In 
the present flight tests, the measured waveforms show 
a relative phase lag that varies from 17° to 47°. The 
minimum lag, 17°, occurs at x/c — 0.4 whereas near 
the leading edge, x/c = 0.05, the lag is about 37° and 
near the trailing edge, x/c = 0.8, it is 47°. 

On the pressure surface (Fig. 2) the predicted 
and measured waveforms differ significantly. This is 
due to the appearance of a computed shock wave dur- 
ing the highly loaded portion of the revolution. The 
measurements do not show such a shock wave. (Such 
a discrepancy was also observed in one of the steady 
(zero angle of attack) Euler computations reported 7 ). 
The measured waveforms at x/c = 0.15 and 0.8 show 
a phase lag compared to the predictions whereas at 
x/c = 0.4 and 0.6 a phase lead is observed. A more 
detailed study is needed to establish if the appearance 
of the passage shock in the predictions is real and not 
an artifact of the numerical scheme employed. 

Figure 3 shows the blade pressure waveforms at 
the transducer locations on the suction surface at the 
outboard radial station, r/R — 0.95. Near the leading 
edge, the peak amplitude is underpredicted. The rel- 
ative phase lag of the measured waveform compared 
to the computed ones exist at all transducer locations 
except that at x/c = 0.58 and it is in the range of 34° 
to 47°. At x/c — 0.58 the measured waveform shows a 
small relative phase lead compared to the predictions. 

The blade pressure waveforms” on - the pressure 
surface at the outboard radial station are shown in 
Fig. 4. At these transducer locations, the magni- 
tudes are substantially overpredicted. Also, signif- 
icant differences in shape of the waveforms are ob- 
served. However, a waveform of the shape similar to 
that predicted at x/c = 0.25 has been observed in 
the wind tunnel tests of Heidelberg 6 on the pressure 
surface for an inflow angle of 2° (for example, see Fig. 
5c in 6 ). However, it has not been possible to estab- 
lish clearly, if the differences between the predicted 
and measured waveforms stem from the installation 
effects or not, due to lack of similar wind tunnel data 
at the high inflow angle considered here. At this out- 
board radial station, r/R = 0.95, the value of the rel- 


ative phase lag of the measured waveform also varies 
widely. 

Compressibility Effects: The effect of Mach num- 
ber on unsteady blade pressure is demonstrated in 
Figs. 5 and 6, in terms of blade surface pressure con- 
tours. Figure 5 shows the blade pressure at four az- 
imuthal positions ($ = 0°, 90°, 180° and 270°) on the 
suction surface for all three Mach numbers. The blade 
pressure variation during a revolution is clearly seen 
to be significant at each Mach number, for the high 
inflow angle, 4.75° considered here. The range of pres- 
sure variation during a revolution is of course maxi- 
mum for M = 0.8. Of particular interest is $ = 90° 
position (in the highly loaded part of the revolution) 
where the area of low pressure region increases with 
increase in Mach number. Figure 6 shows the blade 
pressures at the same four azimuthal positions, as in 
Fig. 5, but on the pressure surface for three Mach 
numbers. It is seen that the cyclic variation of the 
blade pressure during a revolution is nearly the same 
for Mach numbers 0.6 and 0.7. But at Mach number 
= 0.8, a shock appears (in the highly loaded region) 
at $ = 90°. The experimental data, however, do not 
show shock formation as discussed below. 

Next, we present the predicted blade pressure 
waveform variations with Mach number along with 
the flight data. For clarity the predicted and mea- 
sured data are presented side by side to show the 
compressibility effect in each case in all the succeed- 
ing figures. Figure 7 shows the predicted and mea- 
sured waveforms at the transducer locations on the 
suction surface at r/R = 0.68. First of all it is ob- 
served that the measured waveforms show a relative 
phase lag compared to the predicted ones at all the 
three Mach numbers. In general, the effect of increas- 
ing Mach number on the waveform is well predicted. 
Exceptions occur near the leading edge, x/c = 0.05, 
where the maximum amplitudes are overpredicted for 
_ M = DJLand 0.7 and at x/c = 0.8 where a kink ap- 
pears in the waveform for M = 0.8. 

The predicted and measured waveforms on the 
pressure surface at r/R = 0.68 are shown in Fig. 8 
for three Mach numbers. The plots show that when 
the Mach number is increased from 0.6 to 0.7 very 
little change occurs in the waveforms, both in predic- 
tion and experiments. However, the magnitudes are 
overpredicted. For M = 0.8, the predicted waveforms 
show sharp changes in magnitudes indicative of the 
presence of a shock (as shown in Fig. 6). But the ex- 
perimental waveforms show that only near the leading 
edge, x/c = 0.15, a noticeable change in waveform oc- 
curs when the Mach number is increased to 0.8. The 
relative phase lag of the measured waveforms persists 
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at all Mach numbers. 

For the outboard radial station, r/R = 0.95, 
a comparison of predicted and measured waveforms 
is presented in Fig. 9. First, it is observed that 
the measured waveforms are distorted from sinusoidal 
form at all transducer locations for all Mach numbers. 
The distortion reduces with increase in Mach num- 
ber. The tip vortex riding over the transducers for 
(highly loaded) part of the revolution may result in 
the measured shape of the waveform. The predictions 
show linear, sinusoidal behavior of the waveform at 
all transducer locations. The magnitudes are under- 
predicted. However, one trend is correctly predicted, 
that is, with increase in Mach number the response 
decreases at x/c = 0.25 and 0.42. 

Euler solutions are known to predict qualitatively 
the formation of tip vortex 11 in steady (zero angle of 
attack) flow, For angular inflow, the measurements 
seem to indicate that the spatial extent of the tip 
vortex on the blade (suction surface) depends on az- 
imuthal position. During highly loaded part of the 
revolution, the tip vortex may extend over the trans- 
ducer at x/c — 0.08, 0.25, 0.42 and 0.58. An adoptive 
grid which can track correctly the extent of the tip 
vortex and its strength at all azimuthal positions may 
be able to improve the predictions in the outboard re- 
gion. 

Figure 10 shows the pressure waveforms on the 
pressure surface at the outboard radial station, r/R 
= 0.95. The measured waveforms do not show much 
variation in absolute magnitude or shape with Mach 
number. That is true in predictions only for Mach 
numbers 0.6 and 0.7. Even at these Mach numbers the 
magnitudes are overpredicted. At M = 0.8, the pre- 
dictions show significant increase in response, which 
is not observed in experiments. In general at all 
Mach numbers, the maximum response occurs near 
the leading edge and the response reduces gradually 
towards the trailing edge, both in measurement and 
prediction. Significant relative phase lag of the mea- 
sured waveform is observed at all transducer locations 
compared to the predicted ones. 

Concluding Remarks 

The effect of compressibility on unsteady 
blade pressures was studied by solving the three- 
dimensional Euler equations. The unsteadiness is due 
to the operation of the propfan at 4.75° angle of at- 
tack. Three Mach numbers, 0.6, 0.7 and 0.8 were con- 
sidered and the predicted waveforms were compared 
with flight data. 

Comparison of predicted and measured wave- 
forms show that in general, the effect of Mach num- 
ber on unsteady blade pressure waveform is correctly 


predicted. The changes in pressure waveforms are 
minimal when the Mach number is increased from 

0.6 to 0.7 (except on the suction surface at r/R = 
0.95). Increasing the Mach number from 0.7 to 0.8 
produces significant differences in predicted pressure 
levels. However, the predicted appearance of a shock 
at this Mach number is not indicated by the exper- 
iments. The measured waveforms show a relative 
phase lag compared to the predicted ones at most 
of the transducer locations. The phase lag is due to 
the installation effect, that is the installed propfan of 
the flight test compared to the propfan alone config- 
uration of the computation. 

The present numerical procedure is unable to 
reproduce the nonlinear variations of the measured 
waveforms. Perhaps an adoptive grid technique which 
captures accurately the extent and strength of the tip 
vortex may improve the prediction in the outboard re- 
gion. When using the predicted (propfan alone) blade 
pressures to compute near-field acoustics for compari- 
son with flight (installed) data, one has to phase shift 
the predicted pressure waveforms to correspond to in- 
stalled measurement. 
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Table I PTA Follow-on Test Cases 3 : 
Predicted and Measured Total 
Power Coefficients 


Run 

No. 

Mach 

No. 

Blade 

Angle 

Power 

Coef. 

Test 

Euler 

Test 

Euler 

140 

0.8 

56.7 

57.0 

1.493 

1.491 

142 

0.7 

56.6 

57.0 

1.533 

1.549 

144 

0.6 

55.5 

57.0 

1.475 

1.386 


(J) -0 0 


Preplan Inflow Angle - 
Aircraft Angle ot Attack + 
Upwash Angle + 

Proplan Nacelle Tift Angle 



90 


180 ° 



Figure 0. Definition Sketch 


6 



Pressure, kPa , , Pressure, kPa 


0 

-8 
8 

0 k. 


-8 

8 


0 


0 

6 

6 


0 L 


-6 


Experimental Data 

x/c=0.05 Euler Predictions 


x/c=0.15 


x/c=0.25 


x/c=0.4 



x/c=0.6 


x/c=0.8 



90 180 270 360 0 90 180 270 

Azimuth Angle, Deg. Azimuth Angle, Deg. 

Fig.l Pressure Waveforms on the Suction Surface at r/R=0.68, M=0.8 
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Fig.2 Pressure Waveforms on the Pressure Surface at r/R=0.68, M=0.8 
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Fig. 3 Pressure Waveforms on the Suction Surface at r/R=0.95, M=0.8 
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Fig. 4 Pressure Waveforms on the Pressure Surface at r/R=0.95, M=0.8 


4 


8 




9 


Fig. 5 Instantaneous Blade Pressure Distribution on the Suction Surface 
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Fig. 6 Instantaneous Blade Pressure Distribution on the Pressure Surface 
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Fig. 7 Predicted and Measured Pressure Waveforms on the Suction Surface 
at r/R=0.68 
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Euler Predictions 
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Fig. 9 Predicted and Measured Pressure Waveforms on the Suction Surface 
at r/R=0.95 
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variations which are not captured by the present numerical procedure. 
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